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Risk factors for seropositivity to bovine 
coronavirus and bovine respiratory syncytial 


virus in dairy herds 


A. Ohlson, C. Heuer, C. Lockhart, M. Travén, U. Emanuelson, S. Alenius 


A cross-sectional study was conducted to evaluate the association between herd-level 
characteristics, management routines and farm-level antibody status of bovine coronavirus 
(BCV) and bovine respiratory syncytial virus (BRSV) in 257 Swedish dairy herds. The possible 
spatial clustering of positive herds compared with negative herds was also investigated. For 
each herd, a pooled milk sample from five primiparous cows was analysed for the presence of 
antibodies to BCV and BRSV. Herd-level information was obtained by a questionnaire. Logistic 
regression was used to evaluate the association between predictors and antibody status to 
BCV and BRSV. Large herd size, being located in southern Sweden, and not providing boots for 
visitors were found to be associated with being antibody-positive to BCV and BRSV. A short 
distance to the nearest cattle herd was an additional risk factor for BCV. One of the studied 
areas was suitable for spatial analysis. Positive herds were not spatially autocorrelated 

when compared with negative herds as estimated by the K-function regarding both BCV 

and BRSV. This indicates that local factors such as daily visiting milk trucks and wild animals 
were unlikely to be important sources of infection in this area. Moran’s / statistics and semi- 
variogram showed no evidence of spatial autocorrelation in the residuals, indicating that 
remaining unidentified factors are not spatially dependent in the areas under study. 


BOVINE coronavirus (BCV) and bovine respiratory syncytial virus 
(BRSV) are two contagious viruses affecting beef and dairy cattle 
worldwide (Clark 1993, Paton and others 1998, Valarcher and Taylor 
2007). The prevalence of antibodies to BCV and BRSV in bulk tank 
milk was 100 per cent in a nationwide study in England and Wales 
(Paton and others 1998). BCV has tropism for both enteric and res- 
piratory tract epithelium, causing winter dysentery in adult cattle, 
diarrhoea in calves, and various degrees of respiratory tract disease 
(Stair and others 1972, Saif and others 1986, Saif 1990, Alenius and 
others 1991). BRSV replicates in the respiratory tract epithelium and 
can cause respiratory signs, fever and emphysema, and can lead to 
secondary bacterial pneumonia and death (Verhoeff and others 1984, 
Viuff and others 1996). In endemic areas the infections mainly affect 
young animals, whereas in non-endemic areas adult cattle may also 
be affected (Alenius and others 1991, Elvander 1996). Once these 
viruses are introduced into susceptible herds, within-herd transmis- 


sion is generally rapid (Verhoeff and others 1984, Alenius and others 
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1991, Hagglund and others 2007). It has been shown that acquired 
antibodies remain detectable for years, even without reinfection 
(Alenius and others 1991, Elvander 1996), whereas maternal anti- 
bodies are detectable for only a few months. To the authors’ knowl- 
edge, there is no evidence that chronic shedders are involved or that 
BCV and BRSV circulate within a herd for longer periods (Van der 
Poel and others 1993). 

Little is known about the transmission routes of BCV and 
BRSV between cattle herds in Sweden or worldwide. Given the 
contagious nature of these viruses, it is important to identify fac- 
tors that may increase herd-level exposure as well as between-farm 
spread of these viruses in order to target control or prevention 
measures properly. Such measures may also assist with preventing 
other infectious diseases from entering the herd. Previous studies, 
conducted in Norway and Sweden, have evaluated the association 
between herd-level characteristics and BCV and BRSV infections in 
dairy herds. The identified risk factors were similar for both viruses. 
Large herd size was found to be a risk factor compared with small 
herd size (Travén and others 1999, Norstrém and others 2000), as 
was artificial insemination (Al) by farm personnel compared with 
Al by external technicians, conventional compared with organic 
management, and the use of free stalls compared with tie stalls 
(Bidokhti and others 2009). 

To the authors’ knowledge, no other studies have been con- 
ducted regarding management practices as risk factors for the intro- 
duction of BCV and BRSV into dairy herds. Additionally, there 
have been no spatial analyses conducted regarding BCV and BRSV 
infections. Spatial analyses are useful for exploring mechanisms 
of geographical clustering and between-herd spreading of these 
diseases. 

The aim of this study was to evaluate the relation between herd- 
level characteristics and management practices and antibody status of 


BCV and BRSV in Swedish cattle herds from April to May 2007. A 
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secondary objective was to investigate whether there were spatial pat- 
terns in the distribution of the two viruses in the study areas. 


Materials and methods 

Study population 

This was a cross-sectional study of 280 dairy herds within seven 
counties in Sweden (Fig 1). The study was conducted from April 
1 to May 31, 2007. Herds were eligible for inclusion in the study if 
they were members of the local livestock association and enrolled in 
the National Animal Disease Recording System (Emanuelson 1988) 
and the Swedish Official Milk Recording Scheme (Olsson and others 
2001). 

Herd selection varied by county for practical reasons. In Uppland, 
Kalmar and Oland the aim was to include the majority of the existing 
dairy herds within a selected sampling area. Herds were included if the 
farmers were willing to participate as they were routinely visited by 
personnel from the local livestock association. In Halland, Gotland, 
Jamtland and Vasterbotten, a convenience sample of 20 herds as indi- 
cator herds was selected: 10 herds with 30 to 80 cows, and 10 herds 
with more than 80 cows. These herds were also sampled by personnel 
from the local livestock association. Finally, the authors included 16 
large herds (>180 cows) that were considered as being managed by 
progressive farmers (farmers with very good herd management skills 
who were keen to apply new recommendations from veterinarians 
and other experts) according to the local veterinarian. These herds 
were distributed throughout Sweden. None of the herds had known 
antibody status to BCV and BRSV in advance of sampling. All herds 
under study were free from bovine viral diarrhoea virus, bovine leu- 
kaemia virus and bovine herpesvirus type 1 and were not vaccinated 
against BCV or BRSV. The point location of the 257 herds that pro- 
vided data for analysis is shown in Fig 2. 


Data collection 

A pooled milk sample from five primiparous cows that were home- 
bred and had calved since the start of January 2007 was obtained from 
each herd as described by Ohlson and others (2009). The sampling 
was performed by veterinarians and technicians from the local live- 
stock association; 10 ml test tubes containing 1.5 mg of the preserva- 
tive agent bronopol were used to collect the samples, which were not 
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diluted or centrifuged and were stored at -20°C until analysis. The 
milk samples were analysed for the presence of IgG antibodies to 
BCV (Alenius and others 1991) and BRSV (Elvander and others 1995) 
by commercially available indirect ELISAs (SVANOVA Biotech). The 
same batch was used for all analyses. The optical density (OD) at 450 
nm was corrected by subtraction of the negative control antigen OD. 
To adjust between day-to-day variations, the authors calculated the 
per cent positivity (PP) as: (corrected OD/positive control corrected 
OD) x 100. Cut-off was set to PP <20, corresponding to a corrected 
OD of 0.20, which is the cut-off recommended by the manufacturer 
for individual milk samples. Sensitivity is estimated to be 84.6 per 
cent for BCV and 94.6 per cent for BRSV, and specificity 100 per cent 
for both (individual samples). The lower sensitivity of the ELISA for 
BCV antibodies is due to difficulty in detecting weakly positive sam- 
ples (Alenius and others 1991), that is, a risk of false-negative herds. 
Herds that have had an outbreak during the past two years, reflected 
by the primiparous cows sampled for this study, should not be false 
negative, however, because recently infected cows are expected to have 
a high antibody titre. 

A questionnaire (available in Swedish from AO) was administered 
to farmers by veterinarians or technicians from the local livestock 
association during the routine farm visits when milk samples were 
collected. In the event that questionnaires were not administered dur- 
ing sampling, they were sent to the farmers by post. To facilitate spa- 
tial analysis, the X and Y coordinates of the herds were obtained from 
the Swedish Board of Agriculture. 


Spatial analyses of 

point data 

The spatial distribution of BCV- and BRSV-positive farms were 
examined for broad patterns and local dependencies in the distribu- 
tion across the studied areas. The Poisson inhomogeneous K-function 
(Ripley 1977, Diggle and Chetwynd 1991) was used to assess sta- 
tistically whether antibody-positive herds were spatially aggregated 
over antibody-negative herds (ie, local dependencies). The analysis 
was repeated for BCV and BRSV, respectively, in all herds and also 
separately for the herds located within the county of Uppland. Spatial 
aggregation at the national level might not give a true reflection of pos- 
sible clustering because the sampled herds were only a small and non- 
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TABLE 1: Number of herds in levels of categorical variables and means of continuous variables (95 per 
cent confidence interval) of herd characteristics and herd routines as risk factors in 257 Swedish dairy 


herds that were antibody positive /negative to bovine coronavirus (BCV) or bovine respiratory syncytial 
virus (BRSV) in May 2007 


tive if the pooled milk samples 
obtained from primiparous cows 
were positive to BCV and BRSV, 
respectively. Previous studies 


show that BRSV and BCV infec- 


Number of herds ; ; : 
Variable Level BCV-negative  BCV-positive BRSV-negative BRSV-positive All herds tions are effectively spread with- 
sree : in the herd (Verhoeff and others 
Abortions, past six Yes 20 86 18 88 106 1984. Alenius and others 1991 
months No 29 121 22 128 150 eane: ? 
Ifyes,number —-'1.4(1.1t0 1.7) 1.9(1.6to2.1) 1.6(1.1t02.0) 1.8(1.6t02.1) 1.8 (1.6 to 2.0) Hagglund and others 2007). A 
Artificial insemination Technicians 33 114 24 123 147 five-cow pooled milk sample 
ua Farm personnel 16 94 16 94 110 therefore increases the power to 
Biosecurity, importance 1 u 0 u 0 1 discriminate antibody-positive 
(1 Not at all, 5 Very vi 0 1 1 0 1 f tive herd a 
important) 3 5 14 3 16 19 ag a ae ame es 
4 20 90 16 94 110 with a single sample. The predic- 
5 23 102 19 106 125 tor variables included a number 
Biosecurity, do visitors 1 0 1 0) 1 1 of herd-level management and 
respect the 2 1 10 2 9 11 A : A A 
routines of the farm? = Y) 32 4 37 41 ed ee re ne Table 
(1Never, 5 Always) 4 26 127 24 129 153 . Given that none or the contin- 
5 B 36 10 39 49 uous variables was linearly relat- 
Boots provided for No or never used 29 136 21 144 165 ed to the logit of the outcome, 
visitors seldom used 1 13 2 22 4 these were categorised accord- 
Sometimes used 5 23 5 23 28 tiie £6 biolonical plausibilip 
Often used i 24 5 26 31 Tie tO ;020 Oelcds Prately ity OF 
Always used 7 ar J 1 18 equal sized groups. Herd size was 
Coats provided for No or never used 9 55 "1 53 64 divided according to quartiles: up 
visitors Seldom used 2 3 1 4 5 to 38, 39 to 59, 60 to 99, and 100 
Someumes Used e 2 : = 22 or more. Distance to the nearest 
Often used 11 54 8 Syl 65 heap ta 
Always used 19 7 16 74 90 cattle herd in kilometres was cat- 
Commercial breeders Yes 3 7 3 7 10 egorised into four groups: up to 
used for calves No 46 201 37 210 247 0.2 km, over 0.2 to 1 km, over 
Ifyes, quarantine Yes 0 0 0 0 0 1 to 2 km, and 2 km or more. 
ie q t 3 q iD Number of visits during the 
Distance to nearest km 2.4(1.6t03.2) 1.2(1.0t01.4) 2.5(1.4t03.6) 1.2(1.1t0 1.4) 1.4 (1.2 to 1.7) & i 
cattle herd past two weeks by technicians 
Equipment shared Yes 27 119 26 120 146 and veterinarians were grouped 
(past two years) No 21 39 13 97 110 into three categories: one, one 
Herd size Number ofcows 70(53to88)  82(74to91)  76(55to97)  81(72to89)  80(72 to 88) to three, and four or more. The 
Lent or borrowed animals Yes 2 13 1 14 15 two variables ‘providing boots 
past two years No 47 195) 3) 203 242 pn } P ae & 
if yes, quarantine ven 0 1 0 1 1 for visitors’ and ‘providing coats 
No 2 12 1 13 14 for visitors’ were dichotomised 
Organic Management Yes 5 6 3 8 11 as 0 No or yes but not always 
; Ne ae ae oy AV Ze used, and 1 Yes and always used. 
Pasture, animals canreach Yes 8 83 5 86 91 Fat  atieadeaeuard ine binkes 
other herds over fence No 41 123} 35 129 164 pals ie au etce ng Lose 
Purchased animals, past Yes 21 115 17 119 136 curity and visitors respecting the 
two years No 28 93 23 98 121 biosecurity routines of the herd 
If yes, season Pasture 3 16 1 18 19 were scored 1 to 5(5 Very impor- 
instable 2 ae $ us a tant/always respecting) and were 
Both 10 59 il 58 69 both ai hi 
If yes, quarantine Yes 1 9 0 10 10 oth grouped into three catego- 
No 48 106 40 109 126 ries: 1 to 3, 4, and 5. The differ- 
Visitors, past two weeks __ Veterinarians 1.0(0.7t01.4) 1.3(1.1t01.6) 1.0(0.6to1.4) 1.3(1.1to1.5) 1.3 (1.1 to 1.5) ent study areas were allocated 
(numbers) Technicians 2.3(1.7t03.0) 2.5(2.1t02.9) 2.0(1.2to2.7) 2.5(2.2t02.9) 2.4 (2.1 to 2.8) 


random part of the total population; therefore, Uppland was selected 
as the study herds represented approximately 85 per cent of all herds in 
this area. For each virus, separate K-functions were generated for posi- 
tive herds (K,..) and negative herds (K,.)- The difference between the 
two metrics as a function of distance was calculated as D(s)=K _, (s)- 
K...(S): The observed difference can be interpreted as a measure of the 
aggregation of antibody-positive herds over and above that observed 
for the antibody-negative herds at relatively small distances (up to 
100 km). To test the hypothesis that there was no aggregation, the 
location of antibody-positive herds was randomly permuted (n=99 
times) and the observed difference function for each permutation was 
calculated (Chetwynd and Diggle 1998). The upper and lower limits 
of the simulations were subsequently plotted to determine whether 
the observed difference function fell outside the limits of the permuta- 
tion envelope. This would indicate significant spatial aggregation of 
antibody-positive herds relative to antibody-negative herds. 


Risk factor analyses 
The outcomes being evaluated were the presence of antibodies to BCV 
or BRSV in the study herds. Herds were defined as antibody-posi- 


into three regions: south, central, 
and north (Table 2). Collinearity 
between the predictors was 
assessed by calculating the Spearman rank correlation, for which r 
<0.60 was considered as little or no collinearity. 

Separate logistic regression models were used to quantify the 
effects of predictors associated with antibody status (positive v nega- 
tive) for BCV and BRSV, respectively. 


logit(y,) = By ag Bx; +o BX ni 


In this equation, the logit of the observed probability of the i herd 
being antibody-positive, p,, was modelled as a function of m herd-level 
predictor variables. 

The association between the herd-level BCV/BRSV status and 
each of the predictors was first evaluated by chi-squared and Fisher’s 
exact tests. Variables with an a level of less than 0.25 were selected 
for inclusion in the logistic regression model. A backward stepwise 
approach was used for model building, with backward elimination 
of non-significant variables at P>0.05. The presence of confounding 
was assessed by examining the effect of each predictor variable on 
the coefficient of other variables in the model by simultaneously 
adding and removing them into and out of the model and examin- 
ing the change in the coefficients of the remaining model variables. 
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TABLE 2: Prevalence and 95 per cent confidence intervals 
(Cl of herds antibody-positive to bovine coronavirus (BCV) and 


bovine respiratory syncytial virus (BRSV) among primiparous cows 
in 257 Swedish dairy herds for each sampled region 


Region County Number of herds BCV (95% Cl) BRSV (95% Cl) 
South Halland 18 os (NA*) 89 (74 to 100) 
Kalmar 44 91 (82 to 100) 89 (79 to 98) 
Oland 50 96 (90 to rah) 98 (94 to 100) 
Gotland 20 90(76to 100) 90 (76 to 100) 
Central Uppland 69 67 (56 to 78) ie (74 to 92) 
North Jamtland 20 70 (49 to 91) 55 (33 to 77) 
Vasterbotten 20 55 (32 to 77) 75 (55 to 95) 
Specialt 16 81 (61 to 100) 75 (53 to 97) 


“Not applicable, all herds antibody positive 
t Herds considered as progressive according to the local veterinarian, distributed 
throughout Sweden 


If there was a change of 25 per cent or more the predictor variable 
was included in the final model. Interactions between the variables 
included in the model process were also evaluated, but no evidence 
of interactions was found. Statistical significance of the model was 
assessed using the Hosmer-Lemeshow goodness-of-fit test and a like- 
lihood ratio test at a significance level of P<0.05. Plots of Pearson 
residuals versus the predicted values were constructed and evaluated 


for eventual outliers. 


The residuals of the two final models were tested for first- and 
second-order spatial patterns using Moran’s / statistic (Moran 1950) 
and plots of semi-variograms (Isaacs and Srivistava 1989). 


TABLE 3: Chi-squared and Fisher’s exact statistic for the association between herd-level 


antibody status to bovine coronavirus (BCV) and bovine respiratory syncytial virus (BRSV) 
and each of the variables considered for inclusion in a logistic regression model 


Moran’s I statistic was used to determine whether there was 
evidence for dependency (autocorrelation) in the residuals between 
neighbouring herds. First-order neighbours were the farms that shared 
a common border with a farm of interest. Second-order neighbours 
were those of the immediate neighbours of a farm of interest. Moran’s 
! was computed using first- to eighth-order neighbourhood defini- 
tions. The relationship between each pair of herds was defined as 1 
Neighbours and 0 Otherwise. In the presence of positive spatial auto- 
correlation in model residuals the computed Moran’s/ statistic will 
be close to 1 and in its absence its value will be close to zero. The 
statistical significance of the observed Moran’s / statistic was assessed 
using a Monte Carlo permutation approach in which all the residual 
values from the final model were randomly assigned to each herd and 
Moran’s [ calculated on each occasion. The observed Moran’s / statistic 
was then ranked among the simulated values. If the observed statistic 
ranked kth among the 999 simulated values the one-sided significance 
level was k/999. A Moran scatter plot (Anselin 1995) provided a visual 
assessment of spatial autocorrelation among herds. Since the distance 
between farms varied throughout the studied areas, it was the authors’ 
opinion that neighbourhood definitions based on adjacency (rather 
than distance) provided a measure of spatial dependence consistent 
across the entire study area. 

Semi-variograms were produced for the herds located in Uppland 
because the study herds represented approximately 85 per cent of the 
existing dairy herds in the restricted sampling area. This plots the semi- 
variance as a function of distance in kilometres between pairs of herds 
(Isaacs and Srivistava 1989). If herds with more similar residuals were 
closer in space than those with less similar residuals, the semi-variance 
would be expected to increase as a function of distance. This would 
indicate the distances at which the residu- 
als were no longer correlated. Significance 
was assessed using Monte Carlo permuta- 
tion tests. If autocorrelation was present 
in the model residuals, this would indicate 
the presence of local factors involved in the 


5 ie is - 
Exposure variable Description Number of farms Positive farms (%) P dynamics of herd-level antibody status to 
BCV BCV and BRSV, which is not accounted 
Attificial insemination” 0: by external technicians 147 77.6 0.10 for by the variables included in the model. 
1: by farm personnel 110 85.5 ; 
Boots provided for visitorst 0: no or yes but not always used 238 82.4 0.034 : Ee = data ra ae cre ne Sia 
1: yes and always used 18 614 in the R statistical package (R Development 
Distance to nearest 1: <0.2 km 51 94.1 0.002 Core Team 2008) or in Intercooled Stata 
cattle herdt 2: >0.2 to1km 107 83.2 (Stata Statistical Software Release 10.0; 
3:>1to2km 56 71.4 StataCorp) 
4:>2km 36 66.7 . 
Herd sizet 1: $38 cows 60 66.7 0.029 | 
2:39 to 59 cows 62 83.9 Results 
3: 60 to 99 cows 65 86.2 Two hundred and fifty-seven of the 280 
: 4: 2100 cows 66 84.9 herds included provided data for analysis. 
MESES : Hs a ae SH Incomplete questionnaires were obtained 
Pasturet 0: cannot reach other herds 164 75 0.001 from some farms, which resulted in data 
4: can reach other herds 91 912 from 245 and 253 herds being available for 
Region 1: south 140 94.3 <0.001 analysis for the BCV and BRSV models, 
2: oie B 67.0 respectively. Twenty-one of the question- 
ean oy mee naires were done by phone interview. All 
oe 257 farm coordinates were available for the 
Boots provided for visitorst 0: no or yes but not always used 238 86.1 0.011 : 
1: yes and always used 18 61.1 spatial analyses. The overall prevalence of 
Commercial breedingt 0: no 10 70.0 0.19 antibody-positive herds for the period April 
; 1: i ee or eae 247 85.4 to May 2007 was 81 per cent (95 per cent 
Equipment 0: O not st are with other herds 110 88.2 0:19 confidence interval [CI] 76 to 86 per cent) 
1: share with other herds 146 82.2 f 
Herd sizet 1: <38 cows 60 76.7 0.26 or BCV and 84.5 per cent (95 per cent Cl 
2:39 to 59 cows 62 37.1 80 to 89 per cent) for BRSV. Prevalence 
3: 60 to 99 cows 65 89.2 by region with 95 per cent Cl] is shown in 
a = 100/cows 66 83.3 Table 2. There was a significant (P<0.001) 
Pasturet 0: cannot reach other herds 164 78.7 0.001 iti larion ber rhode <cée 
4s can each aiker herds af ante positive correlation between antibody sta 
Region’ Tecan 140 92.9 <0.001 tus to BCV and antibody status to BRSV 
2: central 74 81.1 among the 257 herds included: 18 herds 
ee 3: north 44 63.6 were antibody negative to both BCV and 
Veterinary visits past 1:0 visits 88 78.4 0.13 BRSV: 186 herds were antibody positive to 
two weekst 2:1 to 2 visits 85 89.4 : : 
3. 23 visits 93 855 both BCV and BRSV; 22 herds were posi- 


* chi-squared test 
t Fisher's exact test 
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tive to BCV but negative to BRSV; and 31 
herds were negative to BCV but positive to 
BRSV. 
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TABLE 4: Regression coefficients with standard errors (se), odds ratios (OR) and 95 per 
cent confidence intervals (Cl) of OR as estimated in logistic regression models evaluating 


herd-level risk factors for being antibody-positive to bovine coronavirus (BCV) and bovine 
respiratory syncytial virus (BRSV) in Swedish dairy herds in 2007 


to nearest cattle herd’ from the regression 
model. The result remained unchanged, 
demonstrating that distance had not 
accounted for a second-order effect in the 


Variable Level Coefficient (se) P OR (95% Cl) spatial pattern for BCV. 

Sey ogey Discussion 

Boots for visitors 0: no ref ref 1 Resi t I jneeal sick: 
1: yes -0.27 (0.11) 0.073. (0.76 (0.56-1.03) EBLOT Was SOME aes0clated. wd 

Distance to nearest cattle herd 1:<0.2 km ref ref 1 antibody positivity to BCV and BRSV, 
2:>0.2 to1km -1.25 (0.15) 0.074 0.23 (0.07-1.13) with decreasing OR moving from south 
3:>1 to2km “1.94 (0.08) 0.008 (0.14 (0.03-0.60) to north (Table 4). This agrees with pre- 

ee : oe a (0.06) te au (0.02-0.50) vious nationwide prevalence surveys of 
2:39 to 59 0.84 (1.65) 0.10 2.31 (0.84-6.29) BCV and BRSV (Elvander 1996, Travén 
3: 60 to 99 1.47 (2.67) 0.01 4.36 (1.41-13.44) and others 1999). The geographical differ- 
4:2100 1.14 (1.81) 0.039 3.12 (0.85-6.29) ences might be because of a greater trade 

Region 1:south ref ref 1 of animals between herds and/or higher 
2: central miles (0.12) 0.001 ~—0.19(0.08-0.49) eal dener ia th aera A 
3: north -2.30 (0.05) <0.001 ~—(0.10 (0.37-0.28) BRE Censlly NSE SOUT DSi pats: O 

BRSV model Sweden, but other explanations may also 

Boots for visitors 0:no ref ref 1 be possible. The results from the semi- 

1: yes -0.40 (0.14) 0.006 0.67 (0.51-0.89) variograms show, however, that there 

Herd size (Beste ref ref t was no evidence of spatial dependency 
2:39 to 59 0.75 (0.54) 0.16 2.12 (0.74-6.09) ( d-ord ) 100 k Thi 
3: 60 to. 99 1.14 (0.53) 0.044 ~——-3.13 (1.03-9.52) Secon OnCeE Gy ta SU, Te ices 
4: >100 0.50 (0.47) 0.33 1.65 (0.60-4.55) that the unexplained variance, that is, the 

Region 1: south ref ref 1 model residuals, is not due to factors that 
2: central -0.81 (0.49) 0.094. 0.45 (0.17-1.15) are spatially correlated. 

et Bee Wu) e MOTE OTS ABDS 022) Providing boots for visitors was found 
Herd size Lon ee hon en se ane aioe) to be a protective factor for both BCV 
3 (0) f I I 5 14-6. ‘i P 

3: 60 to. 99 1.14 (0.53) 0.044. ~——«3.13 (1.03-9.52) and BRSV antibody status. This was 
4: 2100 0.50 (0.47) 0.33 1.65 (0.60-4.55) not surprising, as biosecurity is the key 


ref Reference 


Kfunctions for all study herds and herds within the Uppland 
region showed that antibody-positive farms were not spatially aggre- 
gated up to a distance of 100 km, indicating that local spread may not 
be an important factor for the transmission of the two viruses across 
the study herds within the Uppland region. 

The results from the univariable screening for each of the pre- 
dictor variables eligible for inclusion in the logistic regression mod- 
els are shown in Table 3. Seven variables were available for BCV 
and seven for BRSV, with region, herd size and ‘providing boots 
for visitors’ being predictors common for both diseases. For herd 
size, P=0.26 for BRSV in the univariate analysis, but was kept for 
model building because it could be an important confounder. The 
odds ratios (ORs) were similar for the two logistic regression models 
(Table 4). Farms in northern Sweden were much less likely to be 
seropositive than those in the south (OR 0.10 and 0.13 for BCV and 
BRSV, respectively). Independent of region, the OR increased with 
herd size except for herds with more than 100 cows (upper quartile), 
which had slightly lower OR compared with the herds with 66 to 
99 cows (third quartile). Providing boots for visitors was a protec- 
tive factor that reduced the OR for BRSV to 0.67 compared with 
not providing boots, and for BCV the OR was decreased to 0.76 
(P=0.07). Despite P=0.07, the authors chose to keep this biologi- 
cally important variable in the model. The OR for BCV reduced as 
the distance to the nearest cattle herd increased. 

The models showed good fit according to the Hosmer-Lemeshow 
goodness-of-fit test with chi-squared (df 8)=4.79, P=0.78 and chi- 
squared (df 6)=1.77, P=0.94 for BCV and BRSV, respectively. Pseudo 
R? was 23.7 per cent for BCV and 14 per cent for BRSV. Plots of 
Pearson residuals versus the predicted values showed no evidence of 
outliers. 

Moran’s / statistic for the model residuals was -0.7845 (P=0.78) 
and -0.048 (P=0.87) for BCV and BRSV, respectively, indicating 
that there was little or no spatial autocorrelation in the residuals of 
both diseases within the studied areas. The semi-variograms also 
confirmed the lack of spatial dependency in model residuals for the 
Uppland area, indicating a lack of evidence for spatial dependence. 
The authors repeated Moran’s / statistic and semi-variogram using 
residuals from the BCV model after excluding the variable ‘distance 


to preventing contagious diseases. The 
fact that the provision of protective coats 
for visitors was not significantly associ- 
ated with antibody status for any of the 
infections can be explained by the fact that most Swedish veteri- 
narians and technicians bring a change of coat for each herd visit, 
whereas boots are only cleaned. Herd size was also significant for 
both BCV and BRSV, with increasing risk as herd size increased. 
This may be explained by a higher frequency of visitors in larger 
herds than in the very small herds; the slightly lower prevalence in 
herds with more than 100 cows versus herds with 66 to 99 cows 
may simply be due to sampling variation in the high-risk herds, as 
prevalence approached 100 per cent in large herds. Another possible 
explanation may be more thorough implementation of biosecurity 
measures in the very large herds. Previous studies reported either 
positive (Travén and others 1999, Norstrém and others 2000) or no 
(Bidokhti and others 2009) correlations between herd size and BCV 
and BRSV status. 

The similarity of the two final models is biologically plausible 
because BCV and BRSV are both highly contagious diseases and there 
was also a significant positive relationship between antibody status to 
BCV and BRSV. The additional significant factor for BCV, ‘distance 
to nearest cattle herd’, was not significant for BRSV. This could be 
explained by the fact that BCV is shed via faeces, which might be 
more easily spread between herds than nasal discharge, the primary 
means of transmission for BRSV (Clark 1993, Van der Poel and others 
1994). 

The results from analysing point data showed no evidence 
that antibody-positive herds were spatially aggregated over and 
above negative herds, either for BCV or BRSV. This result seems 
to contradict ‘distance to nearest cattle herd’ for BCV. In the logis- 
tic model, however, the variable ‘distance to nearest cattle herd’ is 
the distance to the closest herd, for which the antibody status is 
unknown (it could be a herd not included in the study, for example, 
a beef herd) whereas the K-function is based on the antibody status 
of the neighbours and includes up to the eighth neighbour of the 
study herd. Also, all herds under study were included in the logistic 
model, whereas the spatial analysis of point data included only herds 
from the county of Uppland. The absence of spatial autocorrelation 
regarding both BCV and BRSV in the Uppland area is an important 
finding, showing that local spread may not have a great effect on the 
herd-to-herd transmission of these viruses in areas with moderate 
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animal density. The authors believe that contacts through networks 
(for example, visitors, purchased replacements) could be more impor- 
tant transmission routes for BCV and BRSV in this area. A short 
distance to the nearest cattle herd (less than 200 m) may be related to 
other, unknown, factors. More research is needed, perhaps directed 
at the contact networks between herds, in order to draw further con- 
clusions from this finding. 

A potential weakness of the study is the sampling strategy, 
because the herds included were randomly selected. For Uppland, 
Kalmar and Oland this should not affect the results because the 
majority of the existing herds within limited areas were included. 
In the other counties the selection could lead to a bias, but because 
the antibody status of the herds was unknown before sampling and 
the selection was not based on management skills, the authors do 
not regard this as a major concern. The 16 herds selected to go into 
the special group of large herds, however, could be seriously biased 
by selection, but when running the model without these herds 
the results remained unchanged. There was a risk of information 
bias because farmers had to recall events from the past two years. 
However, most data were collected by interview and this probably 
reduced such a bias to an acceptable level. One strength of this study 
was the high participation rate: 92 per cent of the selected herds 
provided data, and only 5 per cent and 2 per cent of the herds were 
excluded from the final model because of missing values for BCV 
and BRSV, respectively. 

The authors found only one management factor, providing 
boots for visitors, to be associated with antibody status to BCV and 
BRSV. The organic herds did not have a lower prevalence of anti- 
bodies to BCV and BRSV compared with conventional herds, as 
was found to be the case in the study by Bidokhti and others (2009). 
The present study included only 11 organic herds, whereas 20 of the 
40 herds included in the study of Bidokhti and others (2009) were 
organic, and the power of the present study to identify differences 
was consequently lower. The absence of other significant manage- 
ment factors indicates that either non-infected farms were not man- 
aged differently from infected farms, or the questionnaire did not 
include the most critical factors for BRSV and BCV seropositivity. 
That the negative herds managed to stay negative, although being 
surrounded by positive herds, shows that it was possible to avoid 
infection despite the presence of virus in the area. Further inves- 
tigations to elucidate transmission mechanisms, targeting the cen- 
tral and northern parts of Sweden, would be an interesting area for 
future research. 

In conclusion, biosecurity measures in the form of providing boots 
for visitors were associated with a lower herd prevalence of both BCV 
and BRSV in Swedish dairy farms. The prevalence was higher in 
farms located in the southern part of the country than in central or 
northern areas, and increased with herd size. Whereas a short distance 
to the nearest cattle herd increased the probability of seropositivity 
to BCV, spatial analysis, considering the distance between farms as 
well as the infection status of each farm, did not provide any evidence 
that geographical proximity increased the risk of infection with either 
BCV or BRSV. This suggests that local factors such as daily visiting 
milk trucks, wild animals or airborne transmission were unlikely to be 
important sources of infection. 
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